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Abstract. Based on the concept of the cloak generating function, we propose 
a numerical method to compute the invisibility performance of various quadratic 
spherical cloaks and obtain optimal invisibility results. A non- ideal discrete model 
is used to calculate and optimize the total scattering cross-sections of different profiles 
of the generating function. A bell-shaped quadratic spherical cloak is found to be the 
best candidate, which is further optimized by controlling design parameters involved. 
Such improved invisibility is steady even when the model is highly discretized. 
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1. Introduction 

Recently, great progress has been made in botli tlie theory of and experiments on 
invisibility cloaks [H [21 [3], Hj . Wide applications have been found in microwave spectrum 
El El El El [lO] , optical regime [HI [121 [iSl [H [HI HE] , elastodynamics [H [18] , quantum 
mechanics [IHl [20] , and acoustics [211 [221 [231 [21] • One approach to achieve an invisibility 
cloak is to employ the transformation optics (TO) to allow electromagnetic waves 
to be directed around the concealed region and smoothly recovered afterwards. The 
anisotropic parameters of such a cloak are derived from the coordinate transformation. 
This approach was generalized from the cloaking of thermal conductivity [25] and 
then widely applied in many other areas, which provides new approaches to conceal 
passive/active objects [261 [2Z] within their interiors invisible to external illuminations. 
The fundamental idea is the invariance of Maxwell's equations under a space-deforming 
transformation if the material properties are altered accordingly; i.e., a specific spatial 
compression is equivalent to a variation of the material parameters in the flat space. 
Based on TO concept, many efforts have been devoted to the study of 2D cloaks 
(cylindrical [31], elliptical [32], and arbitrary cross-section [33]) due to the simplicity 
in numerical simulations. Inspired by the classic spherical cloak [1], the expressions of 
electromagnetic fields were explicitly presented in terms of spherical Bessel functions 
by Mie theory [26j. However, this analytical scattering theory for classic spherical 
cloaks cannot work if the anisotropic ratio (see the original definition in [3^) is 
anjd;hing other than that in [26]. Two solutions to overcome this problem were 
proposed: 1) multilayers of alternating isotropic layers [35]; and 2) discrete model of the 
inhomogeneous anisotropic shell, and each layer is radially anisotropic but homogeneous 
[36] . Then, the spherical invisibility cloak is near-perfect. Another non-TO route to the 
cloaking in canonical shape is to use a homogenous anisotropic [281 [29] or isotropic 
plasmonic [30] coating. However, in this method, the effectiveness and properties of the 
cloak depend on the object to be cloaked as well as that its size has to be sufficiently 
small compared with the wavelength. Usually, TO-based spherical cloaks need to know 
the prescribed transformation functions first and the required parameters can thus be 
obtained by constructing the explicit transformation matrices. 

Certainly, those sets of parameters from various generating functions for quadratic 
spherical cloaks are ideal, all of which should give zero scattering theoretically. However, 
in actual situations, one has to consider a discrete multilayered model so that the 
invisibility performances of different generating functions distinguish from each other. 
The general method developed in [36] is adopted to calculate the far-field scattering. 
Our numerical results reveal that the power quadratic bell-shaped cloak yields the 
lowest scattering under the same discretization, which is still pronounced when the 
ideal cloaking shell is highly discretized. 

This paper is organized as follows. Section 2 proposes the analysis method which 
will be used to discretize the cloaking shell into multiple layers and then compute the 
far-field diagrams. Section 3 addresses that the bell-shaped profile of the generating 
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function outperforms the other including the hnear one corresponding to the classic 
spherical cloak. Section 4 compares different profiles which give rise to bell-shaped 
profiles. Section 5 discusses the optimization of bell-shaped quadratic cloaks where the 
steady improvement in invisibility performance is verified. 



2. Scattering Algorithm for Arbitrary Spherical Cloaks in Discretized 
Model 

In this section, the scattering theory of multilayer anisotropic spherical particles 
is provided and applied to study a spherical cloak with arbitrary transverse and 
radial parameters. We suppose that the arbitrary field distribution of the incident 
monochromatic wave interacts with the two-layer sphere. The inner sphere is supposed 
to be made of isotropic material. The parameters of the coating depend on the radial 
coordinate and specify the rotationally symmetric anisotropy of the form 

e = er{r)er ^ er + et{r)I, n = iJ,r{r)er ^ er + fit{r)I , (1) 

where Sr and fi,. are the radial dielectric permittivity and magnetic permeability, et and 
fit are the transversal material parameters, / = 1 — e,. Cg) = eg Cg) eg + C?) e^^ is the 
projection operator onto the plane perpendicular to the vector e,., unit vectors e^, e^i, 
and e^p are the basis vectors of the spherical coordinates. 

Using the separation of the variables, the solution of Maxwell's equations in 
spherical coordinates (r, 9, (p) can be presented as 

H{r,e,^)=FUO,^)ii{r), (2) 

where the designation E(r) means that the components of the electric field vector depend 
only on the radial coordinate r as Er{r), Ee{r), and E^{r) (however, the vector itself 
includes the angle dependence in the basis vectors), the second rank tensor in three- 
dimensional space Fim serves to separate the variables {I and m are the integer numbers). 
It can be written as the sum of dyads: 

Flm = Yimer <^ Br + X;^ Bq + (e^ X X;^) (g) e^. (3) 

where Yim{0, ip) and 'K.im{0, ip) are the scalar and vector spherical harmonics. Tensor 
functions are very useful because of their orthogonality conditions 

r r F+^'{0,^)Fim{0,p)smededp = ISmSm'm. (4) 
Jo Jo 

Prom the commutation of the dielectric permittivity £ and magnetic permeability n 
tensors with Fim it follows that the electric and magnetic fields satisfy the system of 
ordinary differential equations 

dH 1 iJlil + 1) 

dr r ^ ^ 



dE , 1 _ V^(^ + l). 



<-r + -<E - ^ ^e^E = ikoptn. (5) 
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where fco = ^^/c the wavenumber in vacuum, uj is the circular frequency of the 
electromagnetic wave. Quantity is called tensor dual to the vector n. It results 
in the vector product, if multiplied by a vector a as n^a = na^ = n x a. 

System ([5]) is the result of the variable separation in Maxwell's equations. Eq. ([5]) 
includes two algebraic scalar equations, therefore, two field components, and E^, can 
be expressed by means of the rest four components. This can be presented as the matrix 
link between the total fields H = Ht + if^Gr and E = Et + E^e^ and their tangential 
components Ht and Et: 

/ 




V{r) 



V 



I 



fj.r(r)kor 



Ht(r) \ 
Et(r) j ' 

\ erir)kor 

Excluding the radial components of the fields from Eq. ([5]), we get to the system of 
ordinary differential equations of the first order for the tangential components joined 
into the four- dimensional vector Wfr) as 



eg 



where 



dW(r) 
dr 



M 



ikoM{r)W{r) 



(7) 




W 




/ He \ 

H^ 



A 



B = Etel 



C = -fite^ + 



/(/ + !) 



Tangential field components Ht and Et play the important part, because they 
are continuous at the spherical interface. Hence, they can be used for solving the 
scattering problem. At first, we will analyze the situation of r-dependent permittivities 
and permeabilities, which arise for the cloak coatings of the spheres. 



Excluding the (^-components of the fields from Eq. 
equation of the second order for wg = e^W = {Hg, Eg 



we derive the differential 







W'g + 



ko^tf^t 





Ml 

fir 



Wg = 0, (9) 



where the prime denotes the r-derivative. Further we will apply the condition on 



the medium parameters, which is usually used for the cloaks: Er 



fir{r) and 



et{r) = /it(r). Then the equations for Hg and Eg coincide and can be written in the 
form 



Wg + 



l{l + 1) St 



ret 



Wg = 0. 



(10) 



This equation can be solved analytically in the very few cases. As an example, we 
can offer the inversely proportional transversal et = ai/r and radial er = a2/r dielectric 
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permittivities. However, such dependencies do not provide the cloak properties of the 
layer. Another analytical solution of Eq. ffTOl) can be obtained for Pendry's cloak, that 
is for the permittivities et = b/{b — a) and Sr = St{r — a)^ jr^ . 

In spite of analytical solutions cannot be found in the all required situations, the 
general structure of solutions can be studied. The solution of two differential equation of 
the second order ([9]) contains four integration constants Ci, C2, d^, and 03. The constants 
can be joined together to a couple of vectors Ci = c\G,q + c'^e,^ and C2 = 026^ + Cge^. 
(^-components of the field vectors and are expressed in terms of the already 
determined ^-components. The link between Q- and (/^-components follows from Eq. 
(I7j). Summing up both components, the resulting field can be presented as 

where ?7i, 7725 Ci^ ^ind C2 are the two-dimensional blocks of the matrix Sir). Quantities 
r^i, Ci) and Ci correspond to the first independent solution of Eq. ([9]), while 172, C25 
and C2 do to the second independent solution. Therefore, the general solution can be 
decomposed into the sum as W = W*^^^ + W^^^ where 

W'"=fM = f'Mc,. w'^)=f H«Uf'Mc. (12) 



Electric and magnetic fields of each independent wave are connected by means of 
impedance tensor F as Etj = Fj-Htj (j = 1,2). Thus, the impedance tensor equals 

r,(r) = 0(r)r/7^(r). (13) 

Vectors Ci and C2 can be expressed by means of the known tangential 

electromagnetic field W(a) as C = S'^^(a)W(a). Then Eq. ffTTl) can be rewritten 
as follows 

W(r) = fi;W(a), = S{r)S~\a), (14) 

where evolution operator (transfer matrix) fi^ connects tangential field components at 
two distinct spatial points, r and a. 

Solution for the fields E(r) and H(r) can be written as the sum over / and m of 
subsequent products of the tensor Fim{0,ip) describing angle dependence (Eq. ([3])), 
matrix V"'(r) restoring the fields with their tangential components (Eq. ([6])), and 
tangential field vectors ffTTl) : 

' Fim{e,^) 

FU9,ip) 

In general, the cloak solutions cannot be studied in the closed form. Therefore, we 
will apply the approximate method of numerical computations, which can be formulated 
as follows. Inhomogeneous spherical shell a < r < b is divided into homogeneous 
spherical layers, i.e. replaced by the multi-layer structure. The number of the layers 
strongly determines the accuracy of calculations, jth homogeneous shell is extended 
from aj^i to aj, where j = 1,...,N, = a and = b. Wave solution of the 
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single homogeneous layer can be presented in the form of evolution operator Qa]-i- The 
solution for the whole inhomogeneous shell is the subsequent product of the elementary 
evolution operators, that is 

The solution of equation with constant permittivities Sr, St and permeabilities 
/ir, /it is expressed by means of the spherical functions. In the layer, the general solution 
is represented using a couple of independent solutions g\f> and g^"*: 

( Heir) \ ^ ( gi'^{hr)c^ + gif{hr)c, \ 

[Ee{r) I [gilKktr)c[ + gif{hr)c'J' ^ ' 

where h = h^euh, = + l)et/ Sr + l/A - 1/2, z/2 = + l)/it//x^ + 1/4 - 1/2. 
Functions g^'"^^ of the order v can be spherical Bessel functions, modified spherical Bessel 
functions, or spherical Hankel functions depending on the problem. 
Blocks 7] and introduced in ( fTTj) are the tensors 

(12) i d(r(/(i'2)) 

Vi,2 = gl^ '^9 ® eg — O e^, 

fit^oT ar 



1 d(rg(^'2)' 
Ci,2 = gil'^^ee ® + ^^"^^^^ ® ^e. 



which can be presented as two-dimensional matrices for computation purposes. 

Now we turn to the scattering of electromagnetic waves from the two-layer sphere. 
The inner sphere of radius a is homogeneous isotropic one {e^^^ and /i*-^-'), while 
the coating a < r < b is characterized by the dielectric permittivity and magnetic 
permeability tensors defined by Eq. ([1]). We suppose that an arbitrary electromagnetic 
field Hinc(r) and Einc(r) is incident onto the two- layer spherical particle from air 

(5(0) = 1, ^(0) = 1). 

Wave solutions in each layer can be written using already known general one (ITSl) . 
Scattered field propagates in air and can be presented in the form of superposition of 
diverging spherical waves. Mathematically, such waves are described by spherical Hankel 
functions of the first kind hl}\x). Denoting the tensors r] and ( with Hankel functions 
gi^^ = h^J^ as f] and (, we get the scattered electromagnetic field 





where F' = C'(^0~^ is the impedance tensor of the Zth scattered wave, Hg™(6) is the 
tangential magnetic field at the particle interface r = b. 

Field inside the isotropic inner sphere is determined by the only spherical Bessel 
function of the first kind. The field in the inhomogeneous shell have no peculiarities and 
can be described in terms of both spherical Bessel functions. Applying the evolution 
operator il^ the electromagnetic field in the shell takes the form 



Hsh(r) \ / Fim \ , , / / 1 
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where r[ = C[{f]{)~^ is the impedance tensor of the Ith wave inside the inner sphere, 
li''{^{a) is the tangential magnetic field at the inner interface of the shell r = a. 

By projecting the fields onto the outer interface r = b and integrating over the angles 



6 and (p with use of orthogonality condition 

/ 



aim 



we derive the boundary conditions 
/ 



where 



W 



Im 



2lT 




(21) 



(22) 



^0 V J'lm[f^^'P)I'^inc{b,0,^) 

Expression (|2T|) represents the system of four linear equations for four components 
of the vectors Hf,™ and H'/". Excluding the constant vector H'/" we derive the amplitude 
of the scattered electromagnetic field 

-1 



t\b) 



T{{a) -/U]^W|™(?3) 



where 



Scattered field can be characterized by the differential cross-section (power radiated 
in e^-direction per solid angle do) 



da 
d^ 



|Hinc(r)p- ^^^^ 
In our notations, the differential cross-section averaged over the azimuthal angle (p 
[over polarizations) takes the form 

2 



da 



sin^d^ 



IH^ 



E 



(25) 



E FiUO,0)Ili^{b) 

l=\m\ 

From the point of view of the scattering theory, it is natural to define the cloak as 
the specially matched layer that provides zero scattering for any material inside. Such 
definition is based on the main property of the cloak — its invisibility (i.e. the cloak 
cannot be detected by optical means). In [] zero scattering was proved analytically for 
the Pendry cloak. 

Zero scattering is specified by the condition Hg™(6) 
using Eq. ( 123|) as follows 



0, which can be rewritten 



T[ia) 



^'6 mc 



0. 



(26) 



Arbitrary incident electromagnetic field W|™ can be excluded from this expression. In 
fact, zero scattering amplitude can be obtained for the trivial situation: electromagnetic 
field is scattered by the air sphere of radius b. This assertion can be presented in the 
form analogous to Eq. ( l26i) : 



( U(6) 



/ W 



Im 



0, 



(27) 



8 



where Fq is the impedance tensor of the Ith wave in the air sphere of radius b. Hence, 
we get to the relation 

(rl(a) -l)nt={T^,ib) -l). (28) 

Impedance tensor T[ contains material parameters e^^'^ and /i*-^-* of the sphere inside 
the cloaking shell. At the same time, zero scattering should be provided for any e'^^^ and 
It can be realized, if the partial derivative of Eq. (!28l) on e^^'^ equals zero, that is 
I ^ = 0. By multiplying this equation by T[ and subtracting it from Eq. ( l28l) . 
we get one more equation, which does not contain the material parameters of the inner 
sphere: (^0 / ^ = — rQ(6) / ^. Finally, we derive the evolution operator of the 
cloaking layer: 

This condition defines the cloak and can be satisfied for the specially chosen 
evolution operator of the cloaking shell. The evolution operator obtained is the 
degenerate block matrix, which inverse matrix is not defined. It should be noted that 
relation (|29l) is independent on the material of the inner sphere. Per se, the derived 
relation connects the wave solutions in the cloak (evolution operator Q'^) and wave 
solutions in the homogeneous air sphere (impedance tensor Fq), that is it performs 
the coordinate transformation for the solutions, but not for the material parameters 
as usually. Unfortunately, it is difficult to determine the dielectric permittivity and 
magnetic permeability of the cloak from Eq. (1291) . 

Cloak condition H'™(6) = substituted to Eq. fl2Tl) results in expression 

fi6W|™= )Hf(a). (30) 

Using (EHD it is clear that H'i™(a) = and E^{^{a) = T{{a)'H.^;^{a) = 0. So, we may 
conclude that both electric and magnetic fields equal zero at the boundary r = a, 
therefore, the electromagnetic field is equal to zero at any spatial point inside the inner 
sphere. In the cloaking shell, the fields equal zero at the inner boundary r = a (owing 
to the continuity of the tangent fields) and equal incident fields at the outer boundary 
r = b. Then, the field inside the cloak is of the form (see Eq. ( 120|) ) 

/ H,h(r) \ ^ f Fim 



Kl(^)^IW|-. (31) 



Esh(r) / fn^'-A Fi 



Im 



3. Bell-Shaped Generating Function for Cloak Optimization 

Starting from this section we consider non-ideal cloaks since we use a discrete model to 
compute and compare the far-field scattering. If the ideal cloaks are considered, each of 
the cloak design is equivalent leading to zero scattering. Realistic cloaks can be made 
of multiple homogeneous spherical layers, which replace the inhomogeneous cloaking 
shell. In this case the scattering is not zero, but noticeably reduced, and such a cloak 
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air 




Figure 1. [Color online] Illustration of the cloaking shell covering the object to be 
concealed. We consider the spherical cloak in free space with the inner radius kQa — tt 
and outer radius knb — 27r. The core material is glass (e'"^"^ — 1.45^ and = 1). 
These quantities are used throughout the whole paper. The material parameters e and 
Ji are determined by applying the proposed transformation-free method to an arbitrary 
cloaking generating function. The cloaking shell is equally divided into A'' layers (each 
layer is homogeneous and anisotropic), and the scattering theory in [36j is used to 
compute the far-field diagrams. 



realization is called non-ideal (see Fig. [T]). In this section we will find the best non-ideal 
cloak providing the lowest cross-section among all designs investigated. 

We will consider some typical generating functions (for transverse dielectric 
permittivities) which exhibit different types of profiles. The simplest generating 
functions are constant, linear, and quadratic ones. Which of them provides the best 
cloaking performance? 

Constant generating function produces the dielectric permittivities of the Pendry's 
classic spherical cloak as it has been demonstrated in the previous section. Linear 
generating function can be generally written as g{r) = r — p, where p is a constant 
parameter. In this case transverse and radial dielectric permittivities become 



Parameter p can take any value except (a + b)/2. It controls the slope of the transverse 
permittivity function. If p < {a + b)/2, et{r) linearly increases, and otherwise it 
monotonically decreases. 

Quadratic generating function has the general form g{r) = {r ~p){r — d) + s, where 
p, d, and s are tunable parameters. The expressions for the permittivities in quadratic 
case can be deduced to 




26(r — p) 



(32) 



b{r — a)'^{r + a — 2p)^ 



(33) 



2r^(r — p){h — a){h + a — 2p) 




b[{r — p){r — d) + s] 



(34) 
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Figure 2. [Color online] Transverse £t and radial dielectric permittivities 
corresponding to different profiles of generating function. Profiles f~9 are described as 
follows: (f ) quadratic generating function with p = 0, d = 6, s = 6^/4; (2) quadratic 
generating function with p = a, d = s ~ {b — a)^/4 + 0.1; (3) quadratic generating 
function with p — a, d = 2b ^ a, s = {b — a)^; (4) linear generating function with 
p ~ a; (5) constant generating function; (6) linear generating function with p = b; (7) 
quadratic generating function with p = a, d = 2b — a, s = 0; (8) quadratic generating 
function with p = a, d = b, s = 0; (9) quadratic generating function with p = 0, d = b, 
s = 0. 



bP'^(r) 

= — , 

r'^[{r — p){r — d) + s]P{b) 

where 



(35) 



P{r) = (p + d) + {pd + s){r-a). (36) 

o ^ 

Quadratic transverse permittivity is a parabola in graphical presentation. The parabola 
can have a minimum (i.e., s > Sq) or maximum (i.e., s < Sq), where Sq = —{b'^ + ab + 
a'^)/3 + {p + d){a + b)/2-pd. 

Using these generating functions, some typical situations depicted in Fig. [2] are 
presented. Profile 5 demonstrates the permittivities for the constant generating function 
g{r) = 1 corresponding to Pendry's cloak. Linear generating functions are presented in 
Profile 4 {g{r) = r — a) and Profile 6 {g{r) = r — b). The other profiles are produced 
using quadratic generating functions. 

The performances of different cloaks can be compared in terms of their scattering 
cross-sections. The best cloaking design possesses the lowest cross-section because of 
the reduced interaction of the electromagnetic wave with the spherical particle. The 
inhomogeneous anisotropic spherical cloaking shell is divided into homogeneous 
anisotropic spherical layers. An experimental realization of this multilayer cloak can 
be the sputtering onto the spherical core. Throughout the whole paper we use = 30. 
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2 4 6 8 10 

Number of profile 

Figure 3. Total cross-sections for profiles of permittivities shown in Fig. [2l The 
number of discrete layers forming the cloak equals = 30. 

In Fig. [3] the total cross-sections resulting from different generating functions are 
shown. Some profiles are approximately equivalent, for example, 1, 2 and 3, or 4 and 6, or 
7 and 9. Profiles 1-3 are characterized by concave- up transverse dielectric permittivity 
e'l > 0. According to Fig. [3] they give rise to the worst results. The flat-curvature 
profiles 4-6 characterized by = are much better. Pendry's cloak (number 5) stands 
out against the other zero-curvature profiles. However, the most effective cloak design is 
the case of concave-down transverse permittivity e" < 0. Profiles 7, 8, and 9 are better 
than profiles 4, 5, and 6, respectively, by approximately 4.8 dB. The quadratic cloak 
with concave-down transverse permittivity (bell-shaped cloak) is shown to be the best 
candidate. 

The maximum of Et in profile 8 is in the middle position of the cloaking shell. 
Shifting the maximum of such a bell shape towards the limit at the outer boundary (i.e., 
profile 7) or inner boundary (i.e., profile 9), the cloaking performance is monotonically 
degraded as shown in Fig. [3l If parameter s is extremely huge in quadratic generating 
function {s oo), the cloak permittivities coincide with those of Pendry's cloak. Thus, 
the increase of s improves the cloak 2 and deteriorates the cloak 8. 

4. The General Class of Bell-Shaped Cloaks 

From the previous section it is concluded that the bell-shaped profile of the transverse 
dielectric permittivity leads to the optimal non-ideal cloaking performance. In the 
present section we will consider the general class of bell-shaped cloaks and choose the 
best type. 

Apart from the quadratic cloak, another three simple bell-shaped profiles will be 
considered: Gaussian, Lorentzian, and Sech. All of them have a single parameter T, 
which sets the width of the profile. We take the maxima of such transverse permittivities 
in the middle of the cloaking shell region (at the point (a -|- b)/2) to compare with 
quadratic cloak. 

Gaussian cloak has the generating function g{r) = exp[— (r — [a + 6)/2)^/(4T^)]. 
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The permittivity functions are 

h 



^* 20FTErf[(6-a)/(4T)]" 

v/^T6(Erf [(6 + a - 2r)/(4T)] - Erf [(6 - a)/(4T)])2^^^^^_ ^^^^ 



2r2Erf [(6 - a)/ {AT)] 

The generating function of the Lorentzian cloak is g{r) = l/[l + {r — {a + b)/2)^/T'^]. 
The transverse and radial permittivities for this cloak are of the form 

b 

^* ~ 2T[1 + (r - (a + b)/2)yT^] arctan[(6 - a)/(2T)] 
_ T6(arctan[(6 + a - 2r)/(2T)] - arctan[(6 - a) /(2T)])2 
~ 2r2 arctan[(6-a)/(2T)] ^ 

(38, 

Sech cloak generating function depends on the radial coordinate as g{r) = sech^[(r — 
(a + b)/2)/T]. The permittivities follows 

6sech2[(r- (a + 6)/2)/T] 



2Ttanh[(6-a)/(2T)] 
T6(tanh[(2r -b- a)/{2T)] - tanh[(6 - a)/{2T)]f 



2r2sech2[(r- (a + 6)/2)/T]tanh[(6-a)/(2T)] ' ^^^^ 
We choose equal 3 dB bandwidths for various transverse permittivity profiles to 
compare different cloaks. Parameters T which are tuned to provide identical 3 dB 
bandwidth for each cloak are given in the caption of Fig. HI In this figure we show 
the total cross-sections of quadratic, Gaussian, Lorentzian, and Sech cloaks. Profiles 
of Gaussian, Lorentzian, and Sech cloaks are very close, resulting in similar scattering 
cross-sections. The influence of the permittivity functions on the cloak performance is 
difficult to tell among these three cloaks. However, it is shown that quadratic cloak in 
Fig. m provides better invisibility when its transverse permittivity vanishes at the inner 
and outer boundaries of the cloaking shell. 

Since the shapes of the Gaussian, Lorentzian, and Sech cloaks are similar, we can 
just select one of them, (e.g., Gaussian) to investigate the significance of the profile, 
which can be varied by the parameter T. The results are demonstrated in Fig. [5l 
The total cross-section has a minimum, which does not provide better cloaking than 
quadratic one though. The cross-section minimization is achieved approximately at 
T = 0.3a. This profile is shown in Fig. [5](a) along with profiles for other T parameters. 
According to this figure the minimization profile has the 3 dB bandwidth equal {b — a)/2. 
Such a profile is neither too narrow nor too wide because narrow profiles {T 0) need 
extremely high discretization and wide profiles (T —>■ oo) tend to the limit of Pendry's 
cloak as shown in Fig. [5](b). 

Thus the bell-shaped quadratic cloak is preferred for non-ideal cloak design, which 
has the lowest cross-section among all bell-shaped cloaks considered in this section. In 
the following section we will show how the quadratic cloak results can be improved. 
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Figure 4. [Color online] (a) Profiles of transverse dielectric permittivity for quadratic 
(profile No. 8 in Fig. Gaussian, Lorentzian, and Sech cloaks and (b) total cross- 
sections of these cloaks. Parameter T equals (fe— a)/4V2 In 2 for Gaussian, (6— a)/(2-\/2) 
for Lorentzian, and (6 — a)/(2-\/21n(\/2 + 1)) for Sech cloak. The number of discrete 
layers forming the cloak equals = 30. 




Figure 5. [Color online] (a) Profiles of transverse dielectric permittivity for Gaussian 
cloaks with different parameters T and (b) total cross-sections versus parameter T. 
The number of discrete layers forming the cloak equals N = 30. 



5. Improved Quadratic Cloaks 

Quadratic cloak is characterized by very simple profile of the transverse dielectric 
permittivity. Also, the quadratic cloak has the scattering almost 5 dB lower than that 
of classic spherical one. Our aim of this section is to find a way of creating the high- 
performance cloaks based on the transformation-free design method and bell-shaped 
quadratic cloak. The high-performance cloak should be similar to the quadratic one. 
Transverse permittivity should have a maximum and vanish at the inner and outer radii 
of the shell: et{a) = etip) = 0. These properties can be satisfied for general generating 
function of the form 



g{r) = (r - a)(r - 6)5-1 (r). 



(40) 
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Figure 6. [Color online] (a) Profiles of transverse permittivity for power quadratic 
cloaks and (b) total cross-sections of these cloaks versus parameter s. Parameters: 
p = a, d = b, N = 30. In (b), only s > is considered for cloaking purposes because 
the total cross-sections corresponding to s < are significantly larger. 



By choosing function gi{r), we can set the permittivity profile of the cloak. The 
function gi{r) can take arbitrary values at the cloak edges r = a and r = b, though it 
should provide the maximum of the transverse permittivity. At first we will consider the 
maximum at the center of the cloak r = {a + b)/2, and then the effect of the non-central 
maximum position will be studied. For instance function gi{r) can be selected with 
Gaussian profile. Then the cloak can be called Gaussian-quadratic one. However, such 
a design is worse than the simple quadratic shape. To provide the better design we will 
focus the quadratic dependence using the gi{r) function 

g,{r) = ((r -p){r-d) + {d- pf/A + s)". (41) 

When n = 0, it is just the bell-shaped quadratic cloak discussed before. The 
permittivities at n > are suppressed due to lengthy expressions. 

In the generating function set by Eqs. fj^Oj) and fHTl) we can vary the power term n 
(the curvature of the transverse permittivity profile at peak), parameters s (the deviation 
from the quadratic cloak) and d (the deviation of the permittivity peak from the center 
of the cloaking region). At n = the generating function is independent on s and d so 
the total cross-section is the straight line in Fig. [6] (the solid line) , where other positive 
values of n are shown as well. The minima of the cross-sections (the best cloaking 
performance) occurs to the parameter s approximately at Smin ~ 9n, i.e., linear to the 
power n. At larger parameter s, the curves tend to the cross-section of the quadratic 
cloak. At small and negative s, the shape of the transverse permittivity contains the 
minimum and a couple of maxima, therefore the total cross-section is substantially 
increased. In Fig. [6](b), the cloaking performance is obviously improved compared with 
the quadratic cloak. Let us further study the effect of the peak position of the profile, 
which is controlled by the parameter d. 

The case d = b describes that the position of the permittivity maximum is in the 
center of the cloaking shell region, li d < b {d > b), the maximum is shifted towards the 
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Radial Coordinate k^r Parameter dib 

Figure 7. [Color online] (a) Profiles of transverse dielectric permittivity for power 
quadratic cloaks and (b) total cross-sections of these cloaks vs. parameter d. 
Parameters: p = a, n = 2, s = 18, N = 30. 

outer (inner) radius of the cloaking shell. Fig. [7] shows that the central position of the 
permittivity maximum is not the optimal choice. The minimization of the cross-section 
is achieved for d ^ 0.846. Such a non-central position is expected to result from the 
spherically curvilinear geometry of the cloak. 

Compared with the quadratic cloak, the improvement of the performance of power 
quadratic cloak is considerable: the total cross-section is further decreased from —54.84 
dB to —57.52 dB. The improvement is caused by the shape of the profile. The profile 
should be parabolic-like with a slightly deformed shape. 

It is also important to consider the differential cross-sections which provides the 
scattering intensity at an arbitrary angle. In Fig. [8] we show the differential cross- 
sections for some typical cloaking designs designed by the transformation-free method 
and considered in a non-ideal situation. The common feature of the cloaks is the reduced 
backscattering. It is seen that the classic spherical cloak is the most visible one. The 
quadratic cloak (blue dotted line) can provide much lower scattering over almost all 
angles compared with Pendry's and Gaussian's bell-shaped cloak. The power quadratic 
cloak is able to further bring down the scattering of the quadratic cloak near the forward 
direction. 

However, one may question that our non-ideal situation may approach to the ideal 
case when the discretization is high (i.e., is much larger than 30). If so, each cloak 
derived from our proposed reversed algorithm should be more and more identical to 
each other. Theoretically, it is true provided that N ^ oo, while the influence of the 
disretization number on the optimization result is still of significant importance in 
practice. We calculate the total cross-sections for different numbers of spherical layers 
in Fig. [9l In general, we observe the conservation of our conclusions on the optimization 
for = 30, except for that the performance of the Gaussian cloak matches with that of 
the quadratic cloak at = 70. When A^ is small, the dependence of scattering reduction 
on the value of A^ is nonlinear. For great number of layers A^ the curves become mostly 
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Figure 8. [Color online] Differential cross-sections of the cloaks derived by the 
proposed transformation-free method. The cloak parameters of each given design have 
been selected to provide the best performance respectively. Parameters: T = 0.3a 
for Gaussian cloak; s = 18, n = 2, p = a, and d = 0.846 for power quadratic cloak; 
N = 30. 
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Figure 9. [Color online] Total cross-sections of different cloak designs versus the 
number of spherical layers N . Parameters of the cloaks are the same as those in Fig. 

m 



linear as shown in Fig. O 
6. Conclusion 

It has been found that the bell-shaped cloaks provide the smallest interaction of the 
cloaking shell with the electromagnetic radiation under the non-ideal situation (i.e., the 
cloaking shell is discretized into layers). Among the bell-shaped cloaks, we have 
compared quadratic, Gaussian, Lorentzian, and Sech cloaks. The last three are very 
similar in profile shape and dependence of controlling parameters. We have concluded 
that the best performance is achieved when the bell-shaped transverse permittivity 
profiles which vanish at the inner and outer radii of the cloaking shell. The simplest 
design of such a type is the quadratic cloak. Improved invisibility performance can be 
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provided by the power quadratic cloak with the maximum of the permittivity profile 
slightly shifted towards the outer boundary. The decrease of cloak's overall scattering 
is about 7.5 dB compared with the classical Pendry's design, and the improvement is 
steady even when the discretization N is quite high. 
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Abstract. Based on the concept of the cloak generating function, we propose 
a numerical method to compute the invisibility performance of various quadratic 
spherical cloaks and obtain optimal invisibility results. A non- ideal discrete model 
is used to calculate and optimize the total scattering cross-sections of different profiles 
of the generating function. A bell-shaped quadratic spherical cloak is found to be the 
best candidate, which is further optimized by controlling design parameters involved. 
Such improved invisibility is steady even when the model is highly discretized. 
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1. Introduction 

Recently, great progress has been made in botli tlie theory of and experiments on 
invisibility cloaks [H [21 [3], Hj . Wide applications have been found in microwave spectrum 
El El El El [lO] , optical regime [HI [121 [iSl [H [HI HE] , elastodynamics [H [18] , quantum 
mechanics [IHl [20] , and acoustics [211 [221 [231 [21] • One approach to achieve an invisibility 
cloak is to employ the transformation optics (TO) to allow electromagnetic waves 
to be directed around the concealed region and smoothly recovered afterwards. The 
anisotropic parameters of such a cloak are derived from the coordinate transformation. 
This approach was generalized from the cloaking of thermal conductivity [25] and 
then widely applied in many other areas, which provides new approaches to conceal 
passive/active objects [261 [2Z] within their interiors invisible to external illuminations. 
The fundamental idea is the invariance of Maxwell's equations under a space-deforming 
transformation if the material properties are altered accordingly; i.e., a specific spatial 
compression is equivalent to a variation of the material parameters in the flat space. 
Based on TO concept, many efforts have been devoted to the study of 2D cloaks 
(cylindrical [31], elliptical [32], and arbitrary cross-section [33]) due to the simplicity 
in numerical simulations. Inspired by the classic spherical cloak [1], the expressions of 
electromagnetic fields were explicitly presented in terms of spherical Bessel functions 
by Mie theory [26j. However, this analytical scattering theory for classic spherical 
cloaks cannot work if the anisotropic ratio (see the original definition in [3^) is 
anjd;hing other than that in [26]. Two solutions to overcome this problem were 
proposed: 1) multilayers of alternating isotropic layers [35]; and 2) discrete model of the 
inhomogeneous anisotropic shell, and each layer is radially anisotropic but homogeneous 
[36] . Then, the spherical invisibility cloak is near-perfect. Another non-TO route to the 
cloaking in canonical shape is to use a homogenous anisotropic [281 [29] or isotropic 
plasmonic [30] coating. However, in this method, the effectiveness and properties of the 
cloak depend on the object to be cloaked as well as that its size has to be sufficiently 
small compared with the wavelength. Usually, TO-based spherical cloaks need to know 
the prescribed transformation functions first and the required parameters can thus be 
obtained by constructing the explicit transformation matrices. 

Certainly, those sets of parameters from various generating functions for quadratic 
spherical cloaks are ideal, all of which should give zero scattering theoretically. However, 
in actual situations, one has to consider a discrete multilayered model so that the 
invisibility performances of different generating functions distinguish from each other. 
The general method developed in [36] is adopted to calculate the far-field scattering. 
Our numerical results reveal that the power quadratic bell-shaped cloak yields the 
lowest scattering under the same discretization, which is still pronounced when the 
ideal cloaking shell is highly discretized. 

This paper is organized as follows. Section 2 proposes the analysis method which 
will be used to discretize the cloaking shell into multiple layers and then compute the 
far-field diagrams. Section 3 addresses that the bell-shaped profile of the generating 
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function outperforms the other including the hnear one corresponding to the classic 
spherical cloak. Section 4 compares different profiles which give rise to bell-shaped 
profiles. Section 5 discusses the optimization of bell-shaped quadratic cloaks where the 
steady improvement in invisibility performance is verified. 



2. Scattering Algorithm for Arbitrary Spherical Cloaks in Discretized 
Model 

In this section, the scattering theory of multilayer anisotropic spherical particles 
is provided and applied to study a spherical cloak with arbitrary transverse and 
radial parameters. We suppose that the arbitrary field distribution of the incident 
monochromatic wave interacts with the two-layer sphere. The inner sphere is supposed 
to be made of isotropic material. The parameters of the coating depend on the radial 
coordinate and specify the rotationally symmetric anisotropy of the form 

e = er{r)er ^ er + et{r)I, n = iJ,r{r)er ^ er + fit{r)I , (1) 

where Sr and fi,. are the radial dielectric permittivity and magnetic permeability, et and 
fit are the transversal material parameters, / = 1 — e,. Cg) = eg Cg) eg + C?) e^^ is the 
projection operator onto the plane perpendicular to the vector e,., unit vectors e^, e^i, 
and e^p are the basis vectors of the spherical coordinates. 

Using the separation of the variables, the solution of Maxwell's equations in 
spherical coordinates (r, 9, (p) can be presented as 

H{r,e,^)=FUO,^)ii{r), (2) 

where the designation E(r) means that the components of the electric field vector depend 
only on the radial coordinate r as Er{r), Ee{r), and E^{r) (however, the vector itself 
includes the angle dependence in the basis vectors), the second rank tensor in three- 
dimensional space Fim serves to separate the variables {I and m are the integer numbers). 
It can be written as the sum of dyads: 

Flm = Yimer <^ Br + X;^ Bq + (e^ X X;^) (g) e^. (3) 

where Yim{0, ip) and 'K.im{0, ip) are the scalar and vector spherical harmonics. Tensor 
functions are very useful because of their orthogonality conditions 

r r F+^'{0,^)Fim{0,p)smededp = ISmSm'm. (4) 
Jo Jo 

Prom the commutation of the dielectric permittivity £ and magnetic permeability n 
tensors with Fim it follows that the electric and magnetic fields satisfy the system of 
ordinary differential equations 

dH 1 iJlil + 1) 

dr r ^ ^ 



dE , 1 _ V^(^ + l). 



<-r + -<E - ^ ^e^E = ikoptn. (5) 
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where fco = ^^/c the wavenumber in vacuum, uj is the circular frequency of the 
electromagnetic wave. Quantity is called tensor dual to the vector n. It results 
in the vector product, if multiplied by a vector a as n^a = na^ = n x a. 

System ([5]) is the result of the variable separation in Maxwell's equations. Eq. ([5]) 
includes two algebraic scalar equations, therefore, two field components, and E^, can 
be expressed by means of the rest four components. This can be presented as the matrix 
link between the total fields H = Ht + if^Gr and E = Et + E^e^ and their tangential 
components Ht and Et: 

/ 




V{r) 



V 



I 



fj.r(r)kor 



Ht(r) \ 
Et(r) j ' 

\ erir)kor 

Excluding the radial components of the fields from Eq. ([5]), we get to the system of 
ordinary differential equations of the first order for the tangential components joined 
into the four- dimensional vector Wfr) as 



eg 



where 



dW(r) 
dr 



M 



ikoM{r)W{r) 



(7) 




W 




/ He \ 

H^ 



A 



B = Etel 



C = -fite^ + 



/(/ + !) 



Tangential field components Ht and Et play the important part, because they 
are continuous at the spherical interface. Hence, they can be used for solving the 
scattering problem. At first, we will analyze the situation of r-dependent permittivities 
and permeabilities, which arise for the cloak coatings of the spheres. 



Excluding the (^-components of the fields from Eq. 
equation of the second order for wg = e^W = {Hg, Eg 



we derive the differential 







W'g + 



ko^tf^t 





Ml 

fir 



Wg = 0, (9) 



where the prime denotes the r-derivative. Further we will apply the condition on 



the medium parameters, which is usually used for the cloaks: Er 



fir{r) and 



et{r) = /it(r). Then the equations for Hg and Eg coincide and can be written in the 
form 



Wg + 



l{l + 1) St 



ret 



Wg = 0. 



(10) 



This equation can be solved analytically in the very few cases. As an example, we 
can offer the inversely proportional transversal et = ai/r and radial er = a2/r dielectric 
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permittivities. However, such dependencies do not provide the cloak properties of the 
layer. Another analytical solution of Eq. ffTOl) can be obtained for Pendry's cloak, that 
is for the permittivities et = b/{b — a) and Sr = St{r — a)^ jr^ . 

In spite of analytical solutions cannot be found in the all required situations, the 
general structure of solutions can be studied. The solution of two differential equation of 
the second order ([9]) contains four integration constants Ci, C2, d^, and 03. The constants 
can be joined together to a couple of vectors Ci = c\G,q + c'^e,^ and C2 = 026^ + Cge^. 
(^-components of the field vectors and are expressed in terms of the already 
determined ^-components. The link between Q- and (/^-components follows from Eq. 
(I7j). Summing up both components, the resulting field can be presented as 

where ?7i, 7725 Ci^ ^ind C2 are the two-dimensional blocks of the matrix Sir). Quantities 
r^i, Ci) and Ci correspond to the first independent solution of Eq. ([9]), while 172, C25 
and C2 do to the second independent solution. Therefore, the general solution can be 
decomposed into the sum as W = W*^^^ + W^^^ where 

W'"=fM = f'Mc,. w'^)=f H«Uf'Mc. (12) 



Electric and magnetic fields of each independent wave are connected by means of 
impedance tensor F as Etj = Fj-Htj (j = 1,2). Thus, the impedance tensor equals 

r,(r) = 0(r)r/7^(r). (13) 

Vectors Ci and C2 can be expressed by means of the known tangential 

electromagnetic field W(a) as C = S'^^(a)W(a). Then Eq. ffTTl) can be rewritten 
as follows 

W(r) = fi;W(a), = S{r)S~\a), (14) 

where evolution operator (transfer matrix) fi^ connects tangential field components at 
two distinct spatial points, r and a. 

Solution for the fields E(r) and H(r) can be written as the sum over / and m of 
subsequent products of the tensor Fim{0,ip) describing angle dependence (Eq. ([3])), 
matrix V"'(r) restoring the fields with their tangential components (Eq. ([6])), and 
tangential field vectors ffTTl) : 

' Fim{e,^) 

FU9,ip) 

In general, the cloak solutions cannot be studied in the closed form. Therefore, we 
will apply the approximate method of numerical computations, which can be formulated 
as follows. Inhomogeneous spherical shell a < r < b is divided into homogeneous 
spherical layers, i.e. replaced by the multi-layer structure. The number of the layers 
strongly determines the accuracy of calculations, jth homogeneous shell is extended 
from aj^i to aj, where j = 1,...,N, = a and = b. Wave solution of the 
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single homogeneous layer can be presented in the form of evolution operator Qa]-i- The 
solution for the whole inhomogeneous shell is the subsequent product of the elementary 
evolution operators, that is 

The solution of equation with constant permittivities Sr, St and permeabilities 
/ir, /it is expressed by means of the spherical functions. In the layer, the general solution 
is represented using a couple of independent solutions g\f> and g^"*: 

( Heir) \ ^ ( gi'^{hr)c^ + gif{hr)c, \ 

[Ee{r) I [gilKktr)c[ + gif{hr)c'J' ^ ' 

where h = h^euh, = + l)et/ Sr + l/A - 1/2, z/2 = + l)/it//x^ + 1/4 - 1/2. 
Functions g^'"^^ of the order v can be spherical Bessel functions, modified spherical Bessel 
functions, or spherical Hankel functions depending on the problem. 
Blocks 7] and introduced in ( fTTj) are the tensors 

(12) i d(r(/(i'2)) 

Vi,2 = gl^ '^9 ® eg — O e^, 

fit^oT ar 



1 d(rg(^'2)' 
Ci,2 = gil'^^ee ® + ^^"^^^^ ® ^e. 



which can be presented as two-dimensional matrices for computation purposes. 

Now we turn to the scattering of electromagnetic waves from the two-layer sphere. 
The inner sphere of radius a is homogeneous isotropic one {e^^^ and /i*-^-'), while 
the coating a < r < b is characterized by the dielectric permittivity and magnetic 
permeability tensors defined by Eq. ([1]). We suppose that an arbitrary electromagnetic 
field Hinc(r) and Einc(r) is incident onto the two- layer spherical particle from air 

(5(0) = 1, ^(0) = 1). 

Wave solutions in each layer can be written using already known general one (ITSl) . 
Scattered field propagates in air and can be presented in the form of superposition of 
diverging spherical waves. Mathematically, such waves are described by spherical Hankel 
functions of the first kind hl}\x). Denoting the tensors r] and ( with Hankel functions 
gi^^ = h^J^ as f] and (, we get the scattered electromagnetic field 





where F' = C'(^0~^ is the impedance tensor of the Zth scattered wave, Hg™(6) is the 
tangential magnetic field at the particle interface r = b. 

Field inside the isotropic inner sphere is determined by the only spherical Bessel 
function of the first kind. The field in the inhomogeneous shell have no peculiarities and 
can be described in terms of both spherical Bessel functions. Applying the evolution 
operator il^ the electromagnetic field in the shell takes the form 



Hsh(r) \ / Fim \ , , / / 1 
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where r[ = C[{f]{)~^ is the impedance tensor of the Ith wave inside the inner sphere, 
li''{^{a) is the tangential magnetic field at the inner interface of the shell r = a. 

By projecting the fields onto the outer interface r = b and integrating over the angles 



6 and (p with use of orthogonality condition 

/ 



aim 



we derive the boundary conditions 
/ 



where 



W 



Im 



2lT 




(21) 



(22) 



^0 V J'lm[f^^'P)I'^inc{b,0,^) 

Expression (|2T|) represents the system of four linear equations for four components 
of the vectors Hf,™ and H'/". Excluding the constant vector H'/" we derive the amplitude 
of the scattered electromagnetic field 

-1 



t\b) 



T{{a) -/U]^W|™(?3) 



where 



Scattered field can be characterized by the differential cross-section (power radiated 
in e^-direction per solid angle do) 



da 
d^ 



|Hinc(r)p- ^^^^ 
In our notations, the differential cross-section averaged over the azimuthal angle (p 
[over polarizations) takes the form 

2 



da 



sin^d^ 



IH^ 



E 



(25) 



E FiUO,0)Ili^{b) 

l=\m\ 

From the point of view of the scattering theory, it is natural to define the cloak as 
the specially matched layer that provides zero scattering for any material inside. Such 
definition is based on the main property of the cloak — its invisibility (i.e. the cloak 
cannot be detected by optical means). In [] zero scattering was proved analytically for 
the Pendry cloak. 

Zero scattering is specified by the condition Hg™(6) 
using Eq. ( 123|) as follows 



0, which can be rewritten 



T[ia) 



^'6 mc 



0. 



(26) 



Arbitrary incident electromagnetic field W|™ can be excluded from this expression. In 
fact, zero scattering amplitude can be obtained for the trivial situation: electromagnetic 
field is scattered by the air sphere of radius b. This assertion can be presented in the 
form analogous to Eq. ( l26i) : 



( U(6) 



/ W 



Im 



0, 



(27) 
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where Fq is the impedance tensor of the Ith wave in the air sphere of radius b. Hence, 
we get to the relation 

(rl(a) -l)nt={T^,ib) -l). (28) 

Impedance tensor T[ contains material parameters e^^'^ and /i*-^-* of the sphere inside 
the cloaking shell. At the same time, zero scattering should be provided for any e'^^^ and 
It can be realized, if the partial derivative of Eq. (!28l) on e^^'^ equals zero, that is 
I ^ = 0. By multiplying this equation by T[ and subtracting it from Eq. ( l28l) . 
we get one more equation, which does not contain the material parameters of the inner 
sphere: (^0 / ^ = — rQ(6) / ^. Finally, we derive the evolution operator of the 
cloaking layer: 

This condition defines the cloak and can be satisfied for the specially chosen 
evolution operator of the cloaking shell. The evolution operator obtained is the 
degenerate block matrix, which inverse matrix is not defined. It should be noted that 
relation (|29l) is independent on the material of the inner sphere. Per se, the derived 
relation connects the wave solutions in the cloak (evolution operator Q'^) and wave 
solutions in the homogeneous air sphere (impedance tensor Fq), that is it performs 
the coordinate transformation for the solutions, but not for the material parameters 
as usually. Unfortunately, it is difficult to determine the dielectric permittivity and 
magnetic permeability of the cloak from Eq. (1291) . 

Cloak condition H'™(6) = substituted to Eq. fl2Tl) results in expression 

fi6W|™= )Hf(a). (30) 

Using (EHD it is clear that H'i™(a) = and E^{^{a) = T{{a)'H.^;^{a) = 0. So, we may 
conclude that both electric and magnetic fields equal zero at the boundary r = a, 
therefore, the electromagnetic field is equal to zero at any spatial point inside the inner 
sphere. In the cloaking shell, the fields equal zero at the inner boundary r = a (owing 
to the continuity of the tangent fields) and equal incident fields at the outer boundary 
r = b. Then, the field inside the cloak is of the form (see Eq. ( 120|) ) 

/ H,h(r) \ ^ f Fim 



Kl(^)^IW|-. (31) 



Esh(r) / fn^'-A Fi 



Im 



3. Bell-Shaped Generating Function for Cloak Optimization 

Starting from this section we consider non-ideal cloaks since we use a discrete model to 
compute and compare the far-field scattering. If the ideal cloaks are considered, each of 
the cloak design is equivalent leading to zero scattering. Realistic cloaks can be made 
of multiple homogeneous spherical layers, which replace the inhomogeneous cloaking 
shell. In this case the scattering is not zero, but noticeably reduced, and such a cloak 



9 



air 




Figure 1. [Color online] Illustration of the cloaking shell covering the object to be 
concealed. We consider the spherical cloak in free space with the inner radius kQa — tt 
and outer radius knb — 27r. The core material is glass (e'"^"^ — 1.45^ and = 1). 
These quantities are used throughout the whole paper. The material parameters e and 
Ji are determined by applying the proposed transformation-free method to an arbitrary 
cloaking generating function. The cloaking shell is equally divided into A'' layers (each 
layer is homogeneous and anisotropic), and the scattering theory in [36j is used to 
compute the far-field diagrams. 



realization is called non-ideal (see Fig. [T]). In this section we will find the best non-ideal 
cloak providing the lowest cross-section among all designs investigated. 

We will consider some typical generating functions (for transverse dielectric 
permittivities) which exhibit different types of profiles. The simplest generating 
functions are constant, linear, and quadratic ones. Which of them provides the best 
cloaking performance? 

Constant generating function produces the dielectric permittivities of the Pendry's 
classic spherical cloak as it has been demonstrated in the previous section. Linear 
generating function can be generally written as g{r) = r — p, where p is a constant 
parameter. In this case transverse and radial dielectric permittivities become 



Parameter p can take any value except (a + b)/2. It controls the slope of the transverse 
permittivity function. If p < {a + b)/2, et{r) linearly increases, and otherwise it 
monotonically decreases. 

Quadratic generating function has the general form g{r) = {r ~p){r — d) + s, where 
p, d, and s are tunable parameters. The expressions for the permittivities in quadratic 
case can be deduced to 




26(r — p) 



(32) 



b{r — a)'^{r + a — 2p)^ 



(33) 



2r^(r — p){h — a){h + a — 2p) 




b[{r — p){r — d) + s] 



(34) 
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Figure 2. [Color online] Transverse £t and radial dielectric permittivities 
corresponding to different profiles of generating function. Profiles f~9 are described as 
follows: (f ) quadratic generating function with p = 0, d = 6, s = 6^/4; (2) quadratic 
generating function with p = a, d = s ~ {b — a)^/4 + 0.1; (3) quadratic generating 
function with p — a, d = 2b ^ a, s = {b — a)^; (4) linear generating function with 
p ~ a; (5) constant generating function; (6) linear generating function with p = b; (7) 
quadratic generating function with p = a, d = 2b — a, s = 0; (8) quadratic generating 
function with p = a, d = b, s = 0; (9) quadratic generating function with p = 0, d = b, 
s = 0. 



bP'^(r) 

= — , 

r'^[{r — p){r — d) + s]P{b) 

where 



(35) 



P{r) = (p + d) + {pd + s){r-a). (36) 

o ^ 

Quadratic transverse permittivity is a parabola in graphical presentation. The parabola 
can have a minimum (i.e., s > Sq) or maximum (i.e., s < Sq), where Sq = —{b'^ + ab + 
a'^)/3 + {p + d){a + b)/2-pd. 

Using these generating functions, some typical situations depicted in Fig. [2] are 
presented. Profile 5 demonstrates the permittivities for the constant generating function 
g{r) = 1 corresponding to Pendry's cloak. Linear generating functions are presented in 
Profile 4 {g{r) = r — a) and Profile 6 {g{r) = r — b). The other profiles are produced 
using quadratic generating functions. 

The performances of different cloaks can be compared in terms of their scattering 
cross-sections. The best cloaking design possesses the lowest cross-section because of 
the reduced interaction of the electromagnetic wave with the spherical particle. The 
inhomogeneous anisotropic spherical cloaking shell is divided into homogeneous 
anisotropic spherical layers. An experimental realization of this multilayer cloak can 
be the sputtering onto the spherical core. Throughout the whole paper we use = 30. 
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Figure 3. Total cross-sections for profiles of permittivities shown in Fig. [2l The 
number of discrete layers forming the cloak equals = 30. 

In Fig. [3] the total cross-sections resulting from different generating functions are 
shown. Some profiles are approximately equivalent, for example, 1, 2 and 3, or 4 and 6, or 
7 and 9. Profiles 1-3 are characterized by concave- up transverse dielectric permittivity 
e'l > 0. According to Fig. [3] they give rise to the worst results. The flat-curvature 
profiles 4-6 characterized by = are much better. Pendry's cloak (number 5) stands 
out against the other zero-curvature profiles. However, the most effective cloak design is 
the case of concave-down transverse permittivity e" < 0. Profiles 7, 8, and 9 are better 
than profiles 4, 5, and 6, respectively, by approximately 4.8 dB. The quadratic cloak 
with concave-down transverse permittivity (bell-shaped cloak) is shown to be the best 
candidate. 

The maximum of Et in profile 8 is in the middle position of the cloaking shell. 
Shifting the maximum of such a bell shape towards the limit at the outer boundary (i.e., 
profile 7) or inner boundary (i.e., profile 9), the cloaking performance is monotonically 
degraded as shown in Fig. [3l If parameter s is extremely huge in quadratic generating 
function {s oo), the cloak permittivities coincide with those of Pendry's cloak. Thus, 
the increase of s improves the cloak 2 and deteriorates the cloak 8. 

4. The General Class of Bell-Shaped Cloaks 

From the previous section it is concluded that the bell-shaped profile of the transverse 
dielectric permittivity leads to the optimal non-ideal cloaking performance. In the 
present section we will consider the general class of bell-shaped cloaks and choose the 
best type. 

Apart from the quadratic cloak, another three simple bell-shaped profiles will be 
considered: Gaussian, Lorentzian, and Sech. All of them have a single parameter T, 
which sets the width of the profile. We take the maxima of such transverse permittivities 
in the middle of the cloaking shell region (at the point (a -|- b)/2) to compare with 
quadratic cloak. 

Gaussian cloak has the generating function g{r) = exp[— (r — [a + 6)/2)^/(4T^)]. 
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The permittivity functions are 

h 



^* 20FTErf[(6-a)/(4T)]" 

v/^T6(Erf [(6 + a - 2r)/(4T)] - Erf [(6 - a)/(4T)])2^^^^^_ ^^^^ 



2r2Erf [(6 - a)/ {AT)] 

The generating function of the Lorentzian cloak is g{r) = l/[l + {r — {a + b)/2)^/T'^]. 
The transverse and radial permittivities for this cloak are of the form 

b 

^* ~ 2T[1 + (r - (a + b)/2)yT^] arctan[(6 - a)/(2T)] 
_ T6(arctan[(6 + a - 2r)/(2T)] - arctan[(6 - a) /(2T)])2 
~ 2r2 arctan[(6-a)/(2T)] ^ 

(38, 

Sech cloak generating function depends on the radial coordinate as g{r) = sech^[(r — 
(a + b)/2)/T]. The permittivities follows 

6sech2[(r- (a + 6)/2)/T] 



2Ttanh[(6-a)/(2T)] 
T6(tanh[(2r -b- a)/{2T)] - tanh[(6 - a)/{2T)]f 



2r2sech2[(r- (a + 6)/2)/T]tanh[(6-a)/(2T)] ' ^^^^ 
We choose equal 3 dB bandwidths for various transverse permittivity profiles to 
compare different cloaks. Parameters T which are tuned to provide identical 3 dB 
bandwidth for each cloak are given in the caption of Fig. HI In this figure we show 
the total cross-sections of quadratic, Gaussian, Lorentzian, and Sech cloaks. Profiles 
of Gaussian, Lorentzian, and Sech cloaks are very close, resulting in similar scattering 
cross-sections. The influence of the permittivity functions on the cloak performance is 
difficult to tell among these three cloaks. However, it is shown that quadratic cloak in 
Fig. m provides better invisibility when its transverse permittivity vanishes at the inner 
and outer boundaries of the cloaking shell. 

Since the shapes of the Gaussian, Lorentzian, and Sech cloaks are similar, we can 
just select one of them, (e.g., Gaussian) to investigate the significance of the profile, 
which can be varied by the parameter T. The results are demonstrated in Fig. [5l 
The total cross-section has a minimum, which does not provide better cloaking than 
quadratic one though. The cross-section minimization is achieved approximately at 
T = 0.3a. This profile is shown in Fig. [5](a) along with profiles for other T parameters. 
According to this figure the minimization profile has the 3 dB bandwidth equal {b — a)/2. 
Such a profile is neither too narrow nor too wide because narrow profiles {T 0) need 
extremely high discretization and wide profiles (T —>■ oo) tend to the limit of Pendry's 
cloak as shown in Fig. [5](b). 

Thus the bell-shaped quadratic cloak is preferred for non-ideal cloak design, which 
has the lowest cross-section among all bell-shaped cloaks considered in this section. In 
the following section we will show how the quadratic cloak results can be improved. 
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Figure 4. [Color online] (a) Profiles of transverse dielectric permittivity for quadratic 
(profile No. 8 in Fig. Gaussian, Lorentzian, and Sech cloaks and (b) total cross- 
sections of these cloaks. Parameter T equals (fe— a)/4V2 In 2 for Gaussian, (6— a)/(2-\/2) 
for Lorentzian, and (6 — a)/(2-\/21n(\/2 + 1)) for Sech cloak. The number of discrete 
layers forming the cloak equals = 30. 




Figure 5. [Color online] (a) Profiles of transverse dielectric permittivity for Gaussian 
cloaks with different parameters T and (b) total cross-sections versus parameter T. 
The number of discrete layers forming the cloak equals N = 30. 



5. Improved Quadratic Cloaks 

Quadratic cloak is characterized by very simple profile of the transverse dielectric 
permittivity. Also, the quadratic cloak has the scattering almost 5 dB lower than that 
of classic spherical one. Our aim of this section is to find a way of creating the high- 
performance cloaks based on the transformation-free design method and bell-shaped 
quadratic cloak. The high-performance cloak should be similar to the quadratic one. 
Transverse permittivity should have a maximum and vanish at the inner and outer radii 
of the shell: et{a) = etip) = 0. These properties can be satisfied for general generating 
function of the form 



g{r) = (r - a)(r - 6)5-1 (r). 



(40) 
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Figure 6. [Color online] (a) Profiles of transverse permittivity for power quadratic 
cloaks and (b) total cross-sections of these cloaks versus parameter s. Parameters: 
p = a, d = b, N = 30. In (b), only s > is considered for cloaking purposes because 
the total cross-sections corresponding to s < are significantly larger. 



By choosing function gi{r), we can set the permittivity profile of the cloak. The 
function gi{r) can take arbitrary values at the cloak edges r = a and r = b, though it 
should provide the maximum of the transverse permittivity. At first we will consider the 
maximum at the center of the cloak r = {a + b)/2, and then the effect of the non-central 
maximum position will be studied. For instance function gi{r) can be selected with 
Gaussian profile. Then the cloak can be called Gaussian-quadratic one. However, such 
a design is worse than the simple quadratic shape. To provide the better design we will 
focus the quadratic dependence using the gi{r) function 

g,{r) = ((r -p){r-d) + {d- pf/A + s)". (41) 

When n = 0, it is just the bell-shaped quadratic cloak discussed before. The 
permittivities at n > are suppressed due to lengthy expressions. 

In the generating function set by Eqs. fj^Oj) and fHTl) we can vary the power term n 
(the curvature of the transverse permittivity profile at peak), parameters s (the deviation 
from the quadratic cloak) and d (the deviation of the permittivity peak from the center 
of the cloaking region). At n = the generating function is independent on s and d so 
the total cross-section is the straight line in Fig. [6] (the solid line) , where other positive 
values of n are shown as well. The minima of the cross-sections (the best cloaking 
performance) occurs to the parameter s approximately at Smin ~ 9n, i.e., linear to the 
power n. At larger parameter s, the curves tend to the cross-section of the quadratic 
cloak. At small and negative s, the shape of the transverse permittivity contains the 
minimum and a couple of maxima, therefore the total cross-section is substantially 
increased. In Fig. [6](b), the cloaking performance is obviously improved compared with 
the quadratic cloak. Let us further study the effect of the peak position of the profile, 
which is controlled by the parameter d. 

The case d = b describes that the position of the permittivity maximum is in the 
center of the cloaking shell region, li d < b {d > b), the maximum is shifted towards the 
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Figure 7. [Color online] (a) Profiles of transverse dielectric permittivity for power 
quadratic cloaks and (b) total cross-sections of these cloaks vs. parameter d. 
Parameters: p = a, n = 2, s = 18, N = 30. 

outer (inner) radius of the cloaking shell. Fig. [7] shows that the central position of the 
permittivity maximum is not the optimal choice. The minimization of the cross-section 
is achieved for d ^ 0.846. Such a non-central position is expected to result from the 
spherically curvilinear geometry of the cloak. 

Compared with the quadratic cloak, the improvement of the performance of power 
quadratic cloak is considerable: the total cross-section is further decreased from —54.84 
dB to —57.52 dB. The improvement is caused by the shape of the profile. The profile 
should be parabolic-like with a slightly deformed shape. 

It is also important to consider the differential cross-sections which provides the 
scattering intensity at an arbitrary angle. In Fig. [8] we show the differential cross- 
sections for some typical cloaking designs designed by the transformation-free method 
and considered in a non-ideal situation. The common feature of the cloaks is the reduced 
backscattering. It is seen that the classic spherical cloak is the most visible one. The 
quadratic cloak (blue dotted line) can provide much lower scattering over almost all 
angles compared with Pendry's and Gaussian's bell-shaped cloak. The power quadratic 
cloak is able to further bring down the scattering of the quadratic cloak near the forward 
direction. 

However, one may question that our non-ideal situation may approach to the ideal 
case when the discretization is high (i.e., is much larger than 30). If so, each cloak 
derived from our proposed reversed algorithm should be more and more identical to 
each other. Theoretically, it is true provided that N ^ oo, while the influence of the 
disretization number on the optimization result is still of significant importance in 
practice. We calculate the total cross-sections for different numbers of spherical layers 
in Fig. [9l In general, we observe the conservation of our conclusions on the optimization 
for = 30, except for that the performance of the Gaussian cloak matches with that of 
the quadratic cloak at = 70. When A^ is small, the dependence of scattering reduction 
on the value of A^ is nonlinear. For great number of layers A^ the curves become mostly 
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Figure 8. [Color online] Differential cross-sections of the cloaks derived by the 
proposed transformation-free method. The cloak parameters of each given design have 
been selected to provide the best performance respectively. Parameters: T = 0.3a 
for Gaussian cloak; s = 18, n = 2, p = a, and d = 0.846 for power quadratic cloak; 
N = 30. 
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Figure 9. [Color online] Total cross-sections of different cloak designs versus the 
number of spherical layers N . Parameters of the cloaks are the same as those in Fig. 

m 



linear as shown in Fig. O 
6. Conclusion 

We have proposed a transformation-free method to obtain required parameters of a 
spherical cloak, based on the concept of the cloaking generating function. It has been 
found that the bell-shaped cloaks provide the smallest interaction of the cloaking shell 
with the electromagnetic radiation under the non-ideal situation (i.e., the cloaking 
shell is discretized into N layers). Among the bell-shaped cloaks, we have compared 
quadratic, Gaussian, Lorentzian, and Sech cloaks. The last three are very similar in 
profile shape and dependence of controlling parameters. We have concluded that the 
best performance is achieved when the bell-shaped transverse permittivity profiles which 
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vanish at the inner and outer radii of the cloaking shell. The simplest design of such 
a type is the quadratic cloak. Improved invisibility performance can be provided by 
the power quadratic cloak with the maximum of the permittivity profile slightly shifted 
towards the outer boundary. The decrease of cloak's overall scattering is about 7.5 dB 
compared with the classical Pendry's design, and the improvement is steady even when 
the discretization N is quite high. 
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